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Hydrodynamic effects in density waves of granular flows 
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: Abstract 

Granular flows in a narrow pipe are studied numerically by the model taking 
account of hydrodynamic effects of fluid surrounding particles. In the simu- 
lations density waves are observed over the wide range of the Stokes number, 
| which represents the inertial effect of particles. The mechanisms of formation 

of the density waves are considered and two types of density waves, which are 
observed in the systems with zero and finite Stokes numbers, are presented. 
Power spectra of density waves obtained from the simulations with non-zero 
. Stokes number show l/f a power law which is similar to that in experiments. 
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Granular materials are not only important in technology and industry, but also are 
very familiar to us. We can see them in daily life, for example, sugar and salt in kitchen, 
sandy beach and soil on which we live. Although these granular materials consist of simple 
components which are granular particles and surrounding fluid, they show a large variety of 
phenomena and attract a lot of interests of scientists During the last decade, many 

people have studied a lot of subjects on granular flows such as vibrating beds fluidized 
beds PHHT2| , hopper flows pipe flows [[H)|-|TP|] , and so on. 



Pipe flow, which we discuss in this paper, seems to be the simplest system of them, 
in which granular particles is falling in a narrow pipe. Even in such a simple system, 
it is observed in experiments that the uniform flow of particles is unstable and density 
waves appear |15|JT6|JT9[| . In addition the fluctuation of density show 1// Q spectra Hl5|Jl9| . 



By numerical simulations assuming mechanical dissipations such as roughness of walls and 
inelasticity of collisions, the creation of density waves |]16|-[T8|] and the appearance of 



spectra [^] were reported. Studies of the clustering behavior of inelastic particles PD|j2l| 
are in the same context. 

However one can expect naturally that the hydrodynamic effects of fluid surrounding 
particles must play an important role in the phenomena. It is sufficient to remember a simple 
example that one particle is falling down by the gravity. Situation is more difficult in the 
pipe flows because the hydrodynamic effects depend on the density of particle significantly 
|2~2fl . Recently the present author and Hayakawa have constructed a model of granular 



particles that takes account of the hydrodynamic effects by applying the method of the 
Stokesian dynamics p3| , and have succeeded to create realistic bubbling and slugging flows 



in fluidized beds [|T^ 

The aim of this paper is to apply this model to the problems of pipe flows and to 
demonstrate the importance of the hydrodynamic effects. 

Since the simulation method used in this paper is described in detail in Ref. | [12]| , we 
briefly summarize the method. Particles are driven by the gravity, hydrodynamic force 
from surrounding fluid and hard-core interaction in direct collisions between particles. The 
equation of motion of particle a is given as 

Stjv^ = F<?> - e 2 , (1) 

where is the velocity of particle a, Ff is the hydrodynamic force, — e z is the unit 
vector directed to the gravity and St = mUo/6TCfia 2 is the Stokes number. In this paper we 
assume that all particles have same radius a and the mass m. In Eq. ([!]), quantities are 
nondimensionalized by the length a and one-particle sedimentation velocity Uq = mg/6irna, 
where g is effective gravitational acceleration corrected by buoyancy and [i is the viscosity of 
the fluid. Direct collisions between particles are considered to be elastic. Particles interact 
with others through the hydrodynamic force which includes the many body effects of particle 
configuration. In low Reynolds number limit, where the viscosity is dominant, F* can be 
written in the matrix form as 

F( Q) = -^R^)-v^, (2) 

P 

where R is so called resistance matrix scaled by 6717m. As described in Ref. |12| , the resistance 
matrix is constructed by the procedure of Stokesian dynamics |23] with periodic boundary 
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condition [24]]. Here we neglect the rotation of particles for simplicity. In addition, we 



approximate (|l|) and @ as follows, 

St>jv^=v[ a) -^ a \ (3) 
-e 2 = £ R M) • vf , (4) 

P 

where v t denotes the terminal velocity. This approximation means that we neglect many 
body effects only in the relaxation process of particle velocity to the terminal velocity. In 
this sense, St' should be treated as the effective Stokes number. 

In the simulation, we use narrow cells with periodic boundary condition. For computa- 
tional efficiency we only discuss the results on monolayer simulations where particles in each 
unit cell can move only on the plane parallel to the direction of gravity and the width of the 
cell is equal to the diameter of particles. In fact there is no qualitative difference between the 
monolayer and full three-dimensional simulations. All simulations presented in this paper 
are calculated as follows. Particles are settled randomly with no initial velocities and the 
evolution obeys the equation of motion and the collision law. In Table | the parameters 
used in the simulations are summarized. 

At first we investigate the effect of St' which is the inertial effect of particles and the only 
parameter in our model. Spatiotemporal patterns of density are shown in Fig. [l] (a),(b) and 
(c), where St' is equal to 0, 10 and 100 respectively. We choose the other parameters as the 
number of particles in the unit cell N = 30 and the cell lengths (L x , L y ,L z ) = (6a, 2a, 80a). 
The figure shows that uniform state is unstable in all case. We can also observed that one 
sharp cluster are formed and the internal motions of particles in the cluster almost freeze in 
St' = and otherwise broader cluster appears. We note that small (large) St' corresponds 
to the fluid with high (low) viscosity like water (air). This means that in not only air but 
also water density waves can be observed. 

Let us consider the mechanisms of dynamics of this behavior. Our model consists of three 
types of mechanisms: inertial effects of particles, hydrodynamics interactions and hard-core 
collisions. From the macroscopic point of view, the inertial effect and the collisions cause 
the advection and the diffusion of particles. Hydrodynamic effect is realized as the mean 
velocity of particles depending on the local density, which is in general decreasing function 
of the density. Therefore only the collisions can have the stabilizing effect on the disturbance 
of density. We define the kinematic regime and the dynamic regime which are correspond 
to the cases where St' is zero and finite respectively. In the kinematic regime (Fig. [I] 
(a)), dynamical effects such as the advection and the collisions vanish and then systems are 
completely governed by (|]). On the other hand in the dynamic regime (Fig. [I] (c)), the 
dynamical effects dominate. The regime in Fig. p] (b) may be the intermediate regime where 
both effects are important. 

We vary parameter N fixing St' = 10 in order to examine the dependence of density 
wave on the particle density. If the number of particles become large (N = 60), internal 
motions of particles is weakened. Thus dynamical effects are suppressed and particles are 
almost frozen in cluster like in the case of St' = 0. If the number of particles becomes 
small (N = 15), which is (d) in Table |, density waves seem to be unstable though they are 
observed. 
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Let us examine the effect of wall confining fluid and particles. For this purpose, we 
introduce the fixed particles which form vertical wall. It is found from the simulation that 
the wall enhances the horizontal motions and the collisional effects when St' is finite. The 
system with wall behaves like that without wall but with larger St'. 

Next we investigate the power spectra obtained from the simulations in Table | (Fig. [|). 
These spectra are calculated as follows. At first we divide horizontally the unit cell into 20 
sub cells and calculate the time series of density in each sub cell. Then we calculate the 
power spectra of each sub cell by a fast Fourier transform with the Parzen window for the 
last 16384 time series of the simulation and average them over 20 sub cells. Frequencies are 
scaled by f which is the lowest frequency. The simulation (e) is presented as the example 
with no density waves for comparison, which correspond to the case of St' = oo and is 
performed without the hydrodynamic effects and the gravity but with the initial velocity as 
their terminal velocities in the case of St' = 10. We can see that power laws are observed 
for the simulations with finite St' in the range between f c = V c /L z and f s = V/4a, where 
V c and V are the mean velocity of the cluster and the particles respectively, L z is system 
size in direction of the gravity and 4a is the length of the sub cells used in the calculation of 
spectra. Therefore f c and f s correspond to the return time of cluster and the characteristic 
time in small scale motion respectively. The exponents in the range are shown in Table O. 

This kind of spectra of density fluctuation have been observed in hourglass twenty years 
ago [15]. In a recent experiment of pipe flow fl9j , similar power spectra have been obtained 



and the high-frequency limit in them may corresponds to f s in this simulation. From the 
meaning of f c and f s , it is suggested that the power law reflects the behavior in relatively 
small scale rather than the cluster itself. The existence of low-frequency limit in Fig. |2| 
may be due to the limit of system size or the periodicity in the simulations. Therefore we 
need to construct the model without periodic boundary condition in order to investigate the 
behavior of larger scale like clusters and the relation to the real systems. 

In conclusion, by means of the numerical simulations |]12| , it is found that in pipe flows the 
hydrodynamic effects of the fluid surrounding particles play an important role and density 
waves appear in the wide range of St. This means that in the water as well as in the air 
density waves can be observed. Two mechanisms to form density waves, which are inertial 
effects of particles and the hydrodynamic interactions, are suggested. In the case of finite 
St', 1/ f a power laws, which may represent the behavior of particles falling from the cluster, 
are observed. 

Finally we note that there are a lot of resemblances between the pipe flows and the traffic 
flows in an expressway. The kinematic waves have been studied as the model of traffic flows 
In addition our model is formally similar to the model of traffic flows presented in Ref. 



and 1/ / power spectra are also observed in the real traffic [p7] and in the numerical 
model f28|. 

The author thanks Hisao Hayakawa and Toshio Tsuzuki for stimulating discussions, 
Yoshi-hiro Taguchi, Akio Nakahara and Mitsugu Matsushita for helpful comments and en- 
couragements and Teruhisa S. Komatsu for useful suggestion on the resemblance to traffic 
flow. Most of simulations in this paper have been done using the facilities of the Supercom- 
puter Center, Institute of Solid State of Physics, University of Tokyo. 
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TABLES 



TABLE I. The parameters of simulations. St' is the effective Stokes number in Eq. (g), N is 
the number of particles in the unit cell, L x ,L y and L z are the lengths of the unit cell. 



St' N L x L y L z 



(a) 0.0 


30 




6 




2 80 


(b) 10.0 


30 




6 




2 80 


(c) 100.0 


30 




6 




2 80 


(d) 10.0 


15 




6 




2 80 


(e) oo 


30 




6 




2 80 


TABLE II. The characteristic velocities and frequencies and the exponents of power laws. V c 


and V is scaled by one-particle sedimentation velocity Uq 


, and f c and f s 


is scaled by the lowest 


frequency /o- 












Vc 


V 


fc 




fs 


a 


(b) 0.12 


0.28 


8 




350 


-1.47 ±3 x 10~ 2 


(c) 0.21 


0.31 


13 




380 


-1.40 ±3 x 10~ 2 


(d) 0.18 


0.43 


11 




510 


-1.32 ±3 x 10~ 2 
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FIGURES 



FIG. 1. Spatiotemporal patterns of density for the first 6000 scaled time with (a) St' = 0,(b) 
St 1 = 10 and (c) St 1 = 100. Other parameters are listed in Table Q. Dark regions correspond to 
high densities. 

FIG. 2. Power spectra S(f) of density fluctuation of the simulations in Table |. These spectra 
are calculated by a fast Fourier transform with the Parzen window. Frequency is scaled by the 
lowest one. Arrows correspond to f c and f s in Table [ij. 
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